WL-TP-92-014 


0 


AD-A259  767 


Transient  Simulation  of  Ram 
Accelerator  Flowfields 


Randall  P.  Drabczuk,  Captain,  USAF 

Wright  Laboratory,  Armament  Directorate 
Analysis  and  Strategic  Defense  Division 
Electromagnetic  Launcher  Technology  Branch 
101  W.  Eglin  Blvd.,  Ste  219 
Eglin  Air  Force  Base  FL  32542-6810 

G.  Rolader 

Science  Applications  Inti.  Corp. 

Shalimar  FL  32547 

S.  Dash,  N.  Sinha,  B.  York 

Science  Applications  Inti.  Corp. 

Ft.  Washington  PA  19034-3211 


JANUARY  1993 

FINAL  PAPER  FOR  PERIOD  JUNE  -  SEPTEMBER  1992 


Approved  for  public  release;  distribution  is  unlimited. 


WRIGHT  LABORATORY,  ARMAMENT  DIRECTORATE 

Air  Force  Materiel  Command  I  United  States  Air  Force  I  Eglin  Air  Force  Base 


98  2  1  071' 


93-01855 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  in  connection  with  a  definitely  Government-related 
procurement,  the  United  States  Government  incurs  no  responsibility  or  any 
obligation  whatsoever.  The  fact  that  the  Government  may  have  formulated  or 
in  any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be 
regarded  by  implication,  or  otherwise  as  in  any  manner  construed,  as  licensing 
the  holder,  or  any  other  person  or  corporation;  or  as  conveying  any  rights  or 
permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in  any 
way  be  related  thereto. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 

The  Public  Affairs  Office  has  reviewed  this  report,  and  it  is  releasable  to  the 
National  Technical  Information  Service  (NTIS),  where  it  will  be  available  to 
the  general  public,  including  foreign  nationals. 

FOR  THE  COMMANDER 


i 


Even  though  this  report  may  contain  special  release  rights  held  by  the 
controlling  office,  please  do  not  request  copies  from  the  Wright  Laboratory, 
Armament  Directorate.  If  you  qu^ify  as  a  recipient,  release  approval  will  be 
obtained  from  the  originating  activity  by  DTIC.  Address  your  request  for 
additional  copies  to: 


Defense  Technical  Information  Center 
Cameron  Station 
Alexandria  VA  22304-6145 


If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing 
list,  or  if  your  organization  no  longer  employs  the  addressee,  please  notify 
WL/MNSH,  Eglin  AFB  FL  32542-6810,  to  help  us  maintain  a  current  mailing 
list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notice 
on  a  specific  document  requires  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  informetion  ts  estimated  to  average  t  f^our  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaimr^  the  data  needed,  and  completing  and  reviewing  the  c^ieaion  of  information.  Ser>d  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
coiiection  of  information,  including  suggestions  for  r^ucing  this  burden,  to  Washington  Headquarters  Services.  Directorate  for  information  Operations  and  Reports.  121S  Jefferson 
Oavis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188).  Washington,  DC  20S03 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPOST  OATE  .  ...  S^REPORT  TYPE  AND  OATES  COVEREO 

January  1993  Final  June  to  September  1992 

4.  TITLE  AND  SUBTITLE 

Transient  Simulation  of  Ram  Accelerator  Flowfields 

5.  FUNDING  NUMBERS 

PEi  61102F 

PRi  2308 

TAi  BW 

WU:  01 

*■  P.  Drabezuk,  Capt,  USAF 

G.  Rolader 

S.  Dash,  N.  Slnha,  B.  York,  SAIC 

7.  gERFOBMINIa  ORGANIZATION  NAME(S)  AND  ^ORESS(ES) 
wright  Laooratory,  Armament  Inrectorate 

Analysis  and  Strategic  Defense  Division 

Electromagnetic  Launcher  Technology  Branch  (WL/MNSH) 

101  W.  Eglin  Blvd.,  Ste  219  together  with  SAIC,  Shalimar 

Eglin  AFB  FL  32S42-6810,  FL  and  Fort  Washington  PA 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AODRESS(ES) 

The  Air  Force  Office  of  Scientific  Research 

Bolling  AFB  DC  20332-6448 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  NUMBER 

WL-TP-92-014 

11.  SUPPLEMENTARY  NOTES 

Approved  by  PA  for  unlimited  release  in  June  1992.  Published  in  the 
AIAA/SAE/ASUE/ASEE  28th  Joint  Propulsion  Conference  proceedings. 

12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

A 

1 13.  ABSTRACT  (Maximum  200  words)  I 

This  paper  describes  the  development  of  an  advanced  computational  fluid  dynamic 
(CFD)  simulation  capability  in  support  of  the  USAF  Armament  Directorate  ram 
accelerator  research  initiative.  The  state-of-the-art  CRAFT  computer  code  has  been 
specialized  for  high  fidelity,  transient  ram  accelerator  simulations  via  inclusion  of  : 
generalized  dynamic  gridding,  solution  adaptive  grid  clustering,  and  high  pressure  ! 
thermo-chemistry.  Selected  ram  accelerator  simulations  are  presented  that  serve  to  ' 
exhibit  the  CRAFT  code  capabilities  and  identify  some  of  the  principle 
research/design  issues.  j 

1 


i 


14.  SUBJECT  TERMS 

Ram  Accelerator 

Hypervelocity  Launcher 

15.  NUMBER  OF  PAGES 

22 

uenerauzed  Dynamic  Gridding  Adaptive  Grid  Clustering 

High  Pressure  Thermo-chemistry 

16.  PRICE  CODE 

17.  SECURITY  aASSIFICATION 
OF  REPORT 

18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

20.  LIMITATION  OF  ABSTRACT 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

SAR 

Standard  form  298  (Rev  2-89) 


t>y  Mo 

?98  10^ 


MSN  7540-01-280-5500 


PREFACE 


This  work  was  funded  by  WL/lfNSH  and  the  Strategic  Defense  Division  of 
the  Wright  Laboratory,  Arwaaent  Directorate  at  Eglin  AFB  FL  under  the  Raa 
Accelerator  Research  Prograa  of  the  Air  Force  Office  of  Scientific  Research. 
Capt  Randall  P.  Drabczuk,  WL/UNSH,  and  personnel  froa  Science  Applications 
International  Corporation  (SAIC)  in  Shaliaar  FL  and  Ft.  Washington  PA 
perforaed  the  work  during  the  period  of  June  to  Septeaber  1992. 

This  technical  paper  was  presented  at  the  29th  JANNAF  Coabustion 
Subcoaaittee  Meeting  in  Haapton  VA  on  19-22  October  1992. 


For 

NTT  3  ilMiii 
WVM'  Ui 
UnKtaownoad 


DTIC  QUALITY  mOPECTED  3 


S'/ - - - 


A'l.'*  I  labll  itr  Codes 
■/▼all  find/or 
'  ■ jp'j',  lal 


iii/iv  (Blank) 


□  □ 
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Wright  Laboratory 
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ABSTRACT 

This  paper  describes  the  development  of  an  advanced  computational  fluid  dynamic  (CFD)  simulaticm  oqnbility 
in  support  of  the  U.S.  Air  Force  Armament  Directorate’s  ram  accelerator  resouch  initiative.  The  state-of-the-art 
CRAFT  computer  code  has  been  specialized  for  high  fidelity,  transient  ram  accelerator  simulations  via  inclusioo 
of  generalized  dynamic  gridding,  solution  adaptive  grid  clustering,  high  pressure  thermochemistry,  etc.  Selected 
ram  accelerator  simulations  are  pre-sented  which  serve  to  exhibit  the  CRAFT  code’s  capabilities  and  identify  some 
of  the  principal  research/design  issues. 


INTRODUCTION 

At  an  earlier  AIAA  meeting,^  a  new  propulsive-oriented  Navier-Stokes  (NS)  code,  CRAFT,  was  introduced 
which  was  an  extension  of  the  TUFF  aerodynamic  NS  code  of  Molvik  and  Merkle.^  Af^licatioos  of  CRAFT  to 
various  combusting,  multiphase  flow  problems  were  exhibited.  In  this  paper,  the  simulation  of  ram  accelerator 
flowfields  with  CRAFT  will  be  discuss^  -  a  role  undertaken  in  support  of  new  research  initiatives  at  the  Armament 
Directorate  of  Wright  Labs,  U.S.  Air  Force. 

The  Armament  Directorate  of  Wright  Laboratories  (WL/MN),  located  at  Eglin  Air  Force  Base,  in  conjunctioo 
with  the  Air  Force  Office  of  Scientific  Research  (AFOSR),  has  initiated  a  multi-year  reseaidi  program  to  advance 
the  understanding  of  the  basic  physics  of  the  ram  accelerator.  The  research  effort  has  three  principal  areas  of  focus: 

(1)  design  and  construction  of  a  large-bore  ram  accelerator  test  focility  at  Eglin  AFB; 

(2)  investigation  of  diagnostic  techniques;  and, 

(3)  develr^ment  and  execution  of  coiiq>utational  fluid  dyruunic  (CFD) 

The  design  of  the  ram  accelerator  research  facilify  to  be  constructed  at  Eglin  AFB  has  been  initiated.  Present  plans 
call  for  the  system  to  mclude  a  conventional  powder  gun  injector,  an  evacuated  vent  chaitiber,  an  dghteen  meter 
long  ram  accelerator  section  with  a  90  mm  bore  diameter,  and  a  catch-tank. 

The  diagnostic  research  will  focus  on  schlieren,  x-ray  shadowgr^rh,  emission  spectroscopy,  and  a  single-pulse 
imaging  (Rayleigh  scattering).  The  schlieren  technique  depends  on  foe  change  in  illumination  due  to  the  changes 
in  refractive  index  (density  gradients),  which  are  made  visible  by  introducing  subsidiary  elements  into  the  optical 
system,  such  as  a  knife  edge,  resulting  in  foe  appearance  of  foe  intensity  variations  in  foe  location  of  the  focussed 
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imige.  The  x>ray  technique  depends  on  the  change  in  the  illumination  due  to  absorption  along  the  line  of  sight  by 
the  gas  constituents,  in  a  manner  similar  to  optical  shadowgraphy.^*^  which  examines  the  intensity  variations  in  a 
plane  which  is  not  the  plane  of  the  focussed  image.  The  interest  in  using  x-rays  derives  from  a  need  to  have  a 
method  that  operates  successfully  at  high  pressures.  X-ray  shadowgnq>hy  has  been  utilized  to  observe  density 
gradients  in  the  flow  as  the  pressure  increases  up  to  values  in  the  kilobar  regime.  Emission  ^wctroscopy  will  be 
employed  to  infer  flowfield  temperature  and  pressure,  as  well  as  coiiq>osition.  An  antici{Mted  diffic^ty  with 
emission  qiectra  is  that  the  opacity  of  the  flow  system  may  increase  rapidly  with  increasing  pressure.  The  rnearch 
will  identify  the  pressure  regime  over  which  this  technique  will  be  useful.  The  fmal  diagnostic  technique  to  be 
investigated  is  single-pulse  plaruu  imaging.  The  primary  focus  will  be  on  Rayleigh  scattering,  but  potentially  both 
Raman  and  Thonqrson  scattering  will  be  utilized. 

High  fidelity  CFD  simulations  will  be  an  integral  part  of  the  overall  ram  accelerator  research  program.  CFD 
simulations  will  be  used  for  conducting  trade-off  studies  to  support  the  in-house  test  fecility  design  ai^  development 
A  second  objective  of  the  computational  fluid  dynamics  ta^  is  to  support  the  development  of  the  diagnostic 
techniques  and  data  mterpretation.  The  final  and  crucial  longer  term  objective  is  to  develop  a  detailed  understanding 
of  the  complex  physical  and  thermochemical  processes  governing  the  ram  accelerator  conc^t  This  paper  will 
provide  a  background  of  the  CFD  technique  being  utilized  to  support  this  research  initiative  and  present  preliminary 
findings. 

Earlier  simulations  of  ram  accelerator  flowfields  have  primarily  emphasized  a  "snapshot*  prediction  of  the  ram 
flowfield,  e.g.,  a  steady-state  simulation  at  a  fixed  trajectory  point.  Such  studies  include  the  work  of  Nusca  at 
BRL,*  Yungster  and  Eberhardt  at  the  University  of  Washington,^***^  Yungster  at  NASA  Lewis,'*’*^  Imlay  et  al.  at 
AMTEC,*^  Hunqihrey  and  Sobota  at  APRI,^^  and  Li,  Kailasanath  and  Oran  at  NRL.*^  The  work  of  Nusca, 
Yungster  et  al.,  and  Imlay  et  al.  has  involved  the  use  of  modem  finite-volume  upwind/implicit  numerics  (e.g., 
higher-order  Roe/TVD).  Nusca  has  used  uncoupled,  global  diemistry  in  an  overiaid  maimer.  Yungster  and  Imlay 
et  al.  have  used  loosely-coupled  (diagonalized)  finite-rate  chemistry.  Li  et  al.  have  utilized  explicit  FCT  numerics 
with  time-split  global  chemistry.  With  regard  to  turbulence  modeling,  Nusca,  Yungster  and  Imlay  et  al.  have 
employed  algebraic  models  (e.g.,  Baldwin-Lomax),  while  Li  et  al.  have  not  utilized  a  turbulence  model  (their  woric 
is  Euler  or  viscous  with  "boosted*  laminar  viscosity). 

A  workshop  was  recently  held  at  AFOSR  with  a  primary  purpose  of  identifying  principal  research  issues  and 
requirements  for  high  fidelity  CFD  simulation  of  ram  accelerator  flowfields.  The  worksht^  group  felt  that  a  full 
transient  simulation  was  requisite  for  studying  the  basic  physics  and  thermochemistry  of  ram  accelerators.  Very 
little  research  has  been  done  in  this  area.  The  only  work  to  date  using  state-of-the-art  upwind/inqilicit  numerics 
to  analyze  transient  bypervelocity  interior  ballistic  flows  is  that  of  Sinba,  Dash  et  al.  as  q>plied  to  Electro-Thenno- 
Chemical  (ETC)  guns.  Alternate  transient  interior  ballistic  flow  simulation  work  for  ram  accelerators,  for  ETC 
gims,  and  for  liquid  propellant  guns,  has  primarily  involved  the  use  of  structured  and  unstructured  explicit  rn^tboda 
(which  cannot  efficiently  treat  thin  boundary  layers).  The  principal  research  areas  identified  at  the  workshop 
included  numerics  (e.g. ,  dynamic  adaptive  gridding  and  advanced  implicit/upwind  algorithms),  turbulence  modeling 
in  a  transient  environment  (e.g.,  LES  sinutlation),  3D  effects  (e.g.,  vortex  shedding  from  fins)  and  conibustion 
phenomenology  in  the  high  Re,  high  pressure  ram  acceleriuor  environment. 

The  methodology  in  CRAFT  is  quite  new  and  reflects  "state-of-the-art*  c^nbilities  in  fliis  country.  Principal 
milestones  in  the  evolution  of  CRAIT  are  summarized  in  Table  1.  Our  work  with  eariier  NS  codes  such  as 
PARCH‘D’**  has  indicated  the  requirement  to  develtq>  specialized  versions  for  specific  applications.*^*^  This 
philosophy  has  been  adhered  to  in  our  work  with  CRAFT.  Spemalized  research  versions  of  CRAFT  developed  to 
date  are  described  in  Table  11.  This  paper  will  review  the  basic  methodology  available  in  CRAFT  (equations  and 
numerics)  supplementing  descriptions  provided  in  earlier  pubiic^ons.  The  problem  specific  features  of  the  RAM 
version  wall  be  emphasized. 


2 


Table  I.  Evolution  of  Methodology  in  CRAFT  Code 


•  UFWIND/IMPLICIT  HNITE-VOLUME  METHODOLOGY  (ROE/TVD)  APPLIED  TO  EXTERNAL  AERODYNAMIC 
FLOWnELDS  —  1986>1988.  Reference:  Welters  end  Thornes,  'Advencee  in  Upwind  Relexetion  Methode,*  Suae  cf  Ae  Art 
Surveys  on  Compuuutonal  Fluid  Meehaniet,  ASME  Pub.,  198S. 

•  EXTENSIONS  OF  ROEOTYD  METHODOLOGY  TO  MULTI-COMPONENT/REAL  GAS  FLOWS  WITH  CHEML 
CAL/TIIERMAL  NONEQUILIBRIUM  ~  198S-1989.  Reference:  Liu  end  Vinokur,  'Upwind  Algoriihme  for  Generel  Tberwo- 
Cbemicel  Nonequilibiium  Flows,*  AIAA  Psper  89-0201,  len.  1989. 

•  EXPLORATORY  STUDIES  OF  ROE/TVD  NUMERICS  FOR  UNSTEADY  FLOWS,  ACCURATE  WAVE  PROPAGATION 
IN  DUCT  (RESONANCE  TUBE)  —  1989.  Reference:  Ridder  end  Beddini,  *Tuiie-Aceunie  Finile-Vobime  Method  for 
Propuleion  Cbember  Flowe,*  AIAA  Peper  89-2554,  July  1989. 

•  EXTENSIONS  OF  ROE/TVD  METHODOLOGY  FOR  PROPULSIVE  FLOWS  WITH  COMBUSTION  CHEMISTRY, 
ADVANCED  TURBULENCE  MODELS  —  1989-1991.  Reference:  Desh,  'Advenced  Computetionel  Models  for  Anelyzinf  Hi(h- 
Specd  Propulsive  Flowfields,*  1990JANNAF  Propulsion  Meeting,  Oct.  1990. 

•  EXTENSIONS  OF  ROE/TVD  METHODOLOGY  FOR  UNSTEADY  RfXTKET  PLUME  PROBLEMS  WITH 
GAS/PARTICLE  NONEQUILIBRIUM  -  1989-1991.  Reference:  Sinbe,  Hoeengedi  end  Desfa,  Tbe  CRAFT  NS  Code  end 
Preliminery  Applicetions  to  Unsteedy,  Reecting,  Multi-Phese  Plume  Flowfield  Problems,*  1991 JANHAF  Plume  Meeting,  Mey  1991. 

•  EXPLORATORY  STUDIES  OF  ROEmO)  METHODOLOGY  FOR  GUN  FLOWHELDS  WITH  DYNAMIC  GRID 
EXTENSIONS  —  1990-1991.  Reference:  Hosengedi  end  Sinbe,  'Development  of  e  Dynemic  Grid  Nevier^tokes  Solver  for 
Computing  Electro-Thermo-Chemieel  Cun  Configurations,*  SAIC/TR-88,  Dec.  1990. 

•  EXTENSIONS  OF  ROEfTVD  METHODOLOGY  FOR  GASH.IQUID  MIXTURES  FOR  ETC  GUN  FLOWFIELD 
SIMULATION  •>  1992.  Reference:  Sinbe,  Hosengedi,  York  end  Desh,  'First  Principlee  Two-Dimensionel  Modeling  of  ETC 
Interior  Bellieiics,'  29A  JANNAF  Combustiort  Meeting,  October  1992. 


Table  II.  Research  Versions  of  CRAFT  NS  Code 


CODE  NAME 

CRAFT/JR 

CRAFT/ETC 

CRAFT/RNP 

CRAFT/LU 

CRAFT/RAM  | 

APPLICATION 

High-Speed  Jet 
Reseereh 

ETC  Cun  Flow- 
fields 

Rocket  Nozzles/ 
Plumes 

Numerical 

Research 

Ram  Accelerator  I 
Flowfields  1 

SPONSOR 

NASA/LeRC, 

HSRP 

BRL 

MICOM 

Internal  Research 

AFOSRft  1 

WL/MNSH  1 

1  EQUATIONS 

|ISI^SS9I 

1D/2D/AXV3D 

ID/2D/AXIOD 

1D/2D/AXI/3D 

1D/2D/AXI/3D  1 

CHEMISTRY 

Perfect  Ges,  Two- 
Streem 

Imperfect  Get, 
Combustion  Chem¬ 
istry,  Multi-Fbaae 
Chemistry 

SPF  Finite-Rale 
Cbemistry,  Two- 
Stream,  Petftet 

Gas 

Perfect  Gas 

Imperfect  Gas, 
Finite-Rale 

Detailed  Kinetics 

TURBULENCE 

kc, 

Compreesibiiity 
Eiocnsions,  ARS 
Extensions 

kc 

ke/Chien, 

Compressibility 

Extensioiis 

ke 

ke/Chien. 

Advanced  Com¬ 
pressibility  Coricc- 
tions  U 

PARTICULATES 

None 

Equilibrated  Mix- 
Urs,  NonequUi- 
brium  Liquid  Wo- 
pcllatxs 

Nonequilibrium 

None 

None  1 

1  SOLUTION 

Inplicit/Uparind 

(Roe/TVD) 

Strongly-Coupled 

Fluid/Spccies/Tbr^ 

b«lencc(8x8 

Block  InverMn), 
Time-Accurste  or 
Time- Asymptotic 

Implicil/Upariiid 
(Roe/TVD)  on  Dy¬ 
namic  Grid,  Time- 
Accurate  Solution, 
Strongly-Coupled 
Equations,  Vari¬ 
able  Matrix  Size 

Implicii/Upsyind 
Gloc/TVD)on 
nxed  Grid,  Time- 
Aeeurtu  or  Tima- 
Aqimplolic, 
Strongly-Coupled 
Equations,  Vari¬ 
able  Matrix  Size 

LU  Upgrade  for 
Robustness,  Fanet 
Convetfence 
(CFL  -  25-50) 

Implicil/Upwind 
(Roe/TVD)  on  Dy¬ 
namic  Grid,  Tuna- 
Accurate  Sohiiion, 
Strongly-Coupled 
Equations,  Vari¬ 
able  Matrix  Size  | 

1  NEW  WORK 

Exploring  New 
Turbulencp 

Models,  New 
Non-Reflective  BC 

IiKlusion  of  Liquid 
Phase,  Gas/Liquid 
Interfoce,  Droplet 
Formation  Model 

Adaptive  Cridding 
for  Unsteady 
Multi-Phase  Flows 

Rewrite  of  Code 
Structure  to  Opti¬ 
mize  LU  Storage 
for  3D 

Adaptive, 

Dynamic  Gridding 
for  Unsteady 

Flows,  LES 
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GAS-PHASE  EQUATIONS 


The  Re3molds-Avenged  Navier-Stokes  equations  integrated  by  the  RAM  version  of  CRAFT  are  listed  in  Table 
m.  They  include  chemiral  species  equations  for  the  species  mass  fractions  p-,,  and  turbulence  model  equations  for 
the  turbulent  kinetic  energy,  k,  and  turbulence  dissipation  rate  parameter,  c.  The  source  terms  in  D  provide  for 
chemical  non-equilibrium  (id,) ,  and  turbulence  non-equilibrium  (E-pc)  where  £  is  the  production  term  for  turbulent 
kinetic  energy.  The  viscous  fluxes  include  full  cross-flow  stress  terms  which  is  requisite  for  solving  3D  problems, 
e.g.,  projectile  with  fins.  A  two-equation  ke  turbulence  modeP^  has  been  incorporated  in  a  stnmgly  coimled, 
inqrlicit  maimer.  Near-wall  effects  have  been  dealt  with  via  the  low  Reynolds  number  formulation  of  ^eon  and 
the  recent  model  of  Rodi.^  High-speed  compressibility  issues  have  be^  accounted  for  via  the  upgrades  of  Dadi 
et  al.^^  based  on  "compressible  dissipation*  concepts  introduced  by  Sarkar^  and  Zeman.^  The  CRAFT  code 
contains  generalized,  multi-component/multi-step,  finite-rate  chemical  kinetics.  The  species  pioduction/destructiim 
rates  are  determined  via  the  law  of  mass  action  and  chemical  reaction  rate  constants  are  expressed  in  Arrhenius 
form.  The  code  accesses  a  JANNAF-based  thermochemical/transport  databank  whose  features  mimic  those  in 
standardized  codes  such  as  SPF.^^  The  databank  contains  thermodynamic  and  transport  data  curve  fits  for  varied 
species  and  rate  constants  for  pertinent  chemical  reactions. 

NUMERICAL  PROCEDURES 

The  CRAFT  code  is  an  extended  propulsive-oriented  derivative  of  the  TUFF  aerodynamic  code  which  was 
formulated  as  a  thin-layer  NS  code  within  a  finite-volume  framework,  with  higher-order  upwind  inviscid  fluxes 
based  on  the  approximate  Riemann  solver  of  Roe.^*  The  numerics  incorporate  TVD  procedures  outlined  by 
Chakravarthy^^  to  eliminate  spurious  oscillations  generally  associated  with  higher  order  upwind  schemes.  The  fluid 
and  species  equations  are  strongly  coupled,  and  the  solution  is  made  fully  implicit  by  linearizing  all  convec- 
tive/difiiisive  fluxes  as  well  as  the  chemical  source  term  (with  respect  to  all  dependent  variables,  not  just  chemical 
species  as  in  other  codes  of  this  type).  An  ADI  splitting  is  enq>loyed.  A  mo^fied  Newton  iteration  is  employed 
to  eliminate  linearization  errors  in  the  iterative  process.  Elimination  of  linearization/factoiizalion  errors  was  recently 
shown  by  Ridder  and  Beddini^^  to  be  essential  for  accurately  representing  complex  wave  propagation  phenomena 
in  a  non-diffiisive  manner,  as  requisite  for  acoustic  applications  and  instability  phenomena. 

The  incorporation  of  nonequilibrium  kinetics  has  followed  the  finite-volume,  fuUy-implicit  philosophy  of  the 
gas  dynamics.  In  departure  from  earlier  generation  explicit  and  loosely  coupled  schemes,  the  species  a^  fluid 
dynamic  equations  are  strongly  coupled  and  simultaneously  solved,  readily  allowing  information  transfer  between 
the  two  sets  of  equations.  In  contrast,  the  earlier  tech^ques  separately  solved  the  two  sets  of  equations  which 
reduces  the  computational  effort  but  severely  inhibits  the  flow  of  information.  This  can  be  a  severe  Handjeap  for 
strongly  combusting  flows,  e.g.,  blast  waves,  gun  diamber  flowfields,  etc.,  uiiere  significant  coupling  exists 
between  the  chemical  release  of  energy  and  flowfield  pressurization.  A  major  breakthrough  has  been  the  pioneering 
work  of  Vinokur  and  Liu^*  in  reformulating  the  original  "perfect  gas"  Roe  scheme  and  the  associated  Reimaan 
problem  to  account  for  multi-component,  real  gas  flow  with  chemical  and  thermal  nonequilibrium. 

The  extension  of  CRAFT  to  include  a  two-equation  ke  turbulence  model  is  a  straightforward  extension  of  the 
baac  fimte-volume  based  upwind/TVD  methodology.  The  turbulence  model  equaticms  are  cast  as  a  part  of  the 
strongly-coupled  equation  set  for  enhanced  accuracy  and  robusmess.  The  implicit  treatment  is  extended  to  the 
turbulence  nonequilibrium  source  terms  to  eliminate  "stiffoess."  This  is  requisite  for  handling  the  low  Reynolds 
number  near  wall  terms.  The  linearization  of  source  terms  is  a  somewhat  inexact  procedure  and  several  strategies 
have  been  developed.  Instead  of  a  brute  force  mathematical  treatment,  the  linearization  is  being  guided  by 
turbulence  behavior  typical  of  low  Reynolds  number  near  wall  regions  and  high  Reynolds  number  five  shear  zones. 
Various  methods  are  being  explored  to  identify  the  optimal  technique  for  propulsive/plume  flowfields.  In  addition, 
the  selection  of  appropriate  background  levels  for  freestream  fluctuation  intensity  and  length  scale  have  been  found 
to  be  sensitive  parameters  and  must  be  judiciously  specified. 
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Table 


CRAFT  NS  Equations 


m. 


1*1/2  “  Ei-l/2)‘*1‘*f 

k*I/2  ■  Gk-,/2)««d’» 
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0 
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0 

'xV*“y^yy*Vy* 

pw 

pOw+IzP 

0 

e 

(e+P)U 

w, 

Pi 

;  E  ■ 

PiO 

;  D  « 

«2 

;  R  - 

p{'xPix*>Viy*“r^i*) 

P2 

PjO 

: 

i 

Pn-l 

p  iC 

"n-l 

P('xPnx*“V»ny**V>M) 

pic 

pCik 

?-p€ 

K'x!£x*“yky*«»JCz) 

P* 

pOc 

^C,P-Cjp*) 

^xx  •  X(u,*Vy*w^*2MU, 
fxy  "  ^yx-  M(«y+v,) 

»^xz  ■  T*x- 

^x  ■  -KTx*»Wix*'^xy*'"yx 

/Jy  -  -KTy*ury*+vryy*wr„ 

/J,  ■  -KT,^UTy,+VT^*WT„ 

X  -  -SV 

IN  THE  ABOVE  EQUATIONS,  ((x,y,4,  «(x,y,z),  r(s.y.z).  REPRESENT  THE  CURVILINEAR  CO(»tDINATES  IN  THE 
TRANSFORMED  COMPUTATIONAL  SPACE,  THE  MEIRICS  t„  t„  t,  ARE  THE  COMPONENTS  OF  THE  CELL  FACE 
NORMAL  7,  POnmNO  IN  THE  POSITIVE  (  DIRECnON  WITH  L&OTH  EQUAL  TO  CELL  FACE  AREA.  SIMILARLY 
THE  OTHER  MEIRIC  QUANTITIES  ARE  THE  COMPONENTS  OF  THE  VECTORS  m  AND  a.  WHICH  ARE  THE  NORMAL 
VECTORS  OF  THE  CORRESPONDINa  m,  AND  a.  THE  FLUID  PROPERTIES  ARE  GIVEN  BY  THE  DENSITY  a. 
CARTESIAN  VELOCITY  COMPONENTS  (ii,  v,  w),  AND  THE  TOTAL  ENERGY  PER  UNIT  VOLUME.  •,  AND  SPECIES 
CONCENTRATIONS,  p;.  THE  TRANSPORT  PROPERTIES  ARE  GIVEN  BY  m  BEING  THE  LAMINAR  VISCOSITY,  y  BEING 
THE  RATIO  OF  SPECIFIC  HEATS,  AND  K  BEING  THE  THERMAL  CONDUCTIVTTY.  THE  TURBULENCE  EQUATIONS 
ARE  DERIVED  IN  TERMS  OF  THE  TURBULENT  KINETIC  ENERGY,  k.  /Un>  THE  DISSIPATION  RATE,  (.  THE 
CONTRAVARIANT  VELOCITY  COMPONENTS  ARE  DEFINED  AS  FOLLOWS; 

0  -  T,U  ♦  TyV  ♦  f,W 
^  •  m,u  ♦  niyV  ♦  ni,w 

>V  -  n,u  ♦  DyV  ♦  n^w 


FOR  TURBULENT  FLOW  SIMULATION,  THE  LAMIN/Ut  VISCOSITY  ASD  THE  THERMAL  CONDUCnVITY  ARE  REPLACED 
BY  THE  RESPECTIVE  SUMS  OF  LAMINAR  AND  TURBULENT  VALUES.  THE  k-<  TURBULENCE  MODEL  CONSTANTS 
WHICH  HAVE  BEEN  USED  ARE  C,  -  .09.  C,  ■  1 .44  «ad  C^  -  1.92. _ 
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DYNAMIC  GRID  UPGRADES 


As  a  preliminary  st^  towards  developing  a  gun  capability,  CRAFT  was  upgraded  to  allow  for  unsteady 
simulations  with  a  moving  projectile  and  thus  an  ever  expanding  and/or  moving  gun  chamber  computational  domain. 
The  generalized  dynamic  grid  formulation  incorporated  enhances  the  capabilities  of  the  code  to  simulate  a  wide 
variety  of  other  unsteady  phenomena,  e.g.,  the  unsteady  interactions  between  deformable  or  non-stadonary 
boundaries  and  fluid  phenomena  which  result  in  time-varying  flowflelds.  Dynamic  grids  are  also  applicable  to 
problems  which  require  a  grid  that  continuously  adapts  to  moving  discontinuities.  This  feature  is  requirite  in 
transient  problems  s^ere  strong  shocks  can  convect  rapidly  through  the  flowfleld. 

The  dynamic  grid  capability  has  been  formulated  in  a  manner  consistent  with  the  finite-volume  philosophy  of 
the  CRAFT  code.  We  have  closely  followed  the  generalized  procedure  outlined  by  Vinokur.^^  Table  IV 
summarizes  this  methodology,  while  additional  details  are  given  in  Ref.  33.  To  briefly  describe  the  procedure,  the 
unsteady  movement  of  the  physical  grid  is  treated  purely  as  a  geometric  quantity.  Therefore,  unlike  a 
finite-difference  formulation  vdiere  the  contravariant  velocity  is  altered  such  that  the  computational  grid  remains 
invariant  in  space  and  time,  the  computational  grid  in  a  finite-volume  formulation  varies  with  time.  The  additional 
flux  due  to  the  motion  of  the  cell  face  is  included  in  the  flux  evaulation  in  an  upwind  fashion.  The  conservati<m 
equations  incorporate  the  change  in  the  flux  resulting  from  the  time-varying  cell  volume  as  a  source  term.  The 
variation  of  the  grid  with  time  has  to  be  specified  depending  on  the  application.  For  a  moving  boundary,  the 
movement  of  the  boundary  has  to  be  related  to  the  force  balance  on  it.  On  the  other  hand,  for  a  grid  adapting  to 
convecting  discontinuities,  the  criteria  for  the  adaption  (e.g.  pressure  or  density  gradient)  must  be  specified. 

RAM  ACCELERATOR  STUDIES 


Steady-State  Simulations 

It  was  imperative  to  first  establish  a  baseline  capability  at  the  steady-state  axisymmetric  level  to  establish  that 
CRAFT  provided  predictions  in  accord  with  those  of  oAer  investigators.*’* With  this  objective  in  mind, 
predictions  were  generated  for  a  projectile  that  was  recently  tested  at  the  University  of  Washington’s  38  mm  ram 
accelerator  test  facility. 

The  projectile  dimensions  are  shown  in  Figure  1.  The  tube  was  pressurized  to  39  atm  with  a  premixed 
combustible  gas  mixture  (2. 75  CH4  -t-  5.8  N2  +  2  Oj).  CRAFT  predictions  were  obtained,  at  the  Euler  non- 
reacting  limit,  for  a  trajectory  point  corresponding  to  a  projectile  flight  Mach  number  of  3.6.  A  301  x  71 
curvilinear  body  fitted  grid  (shown  in  Figure  2)  was  utilized.  Figures  3  and  4  show  computed  pressure  and  Mach 
number  contours,  respectively.  The  overall  flow  structure  is  well  represented  and  in  accord  with  predictions 
obtained  by  the  previous  investigators  utilizing  comparable  numerics.  The  oblique  shock  structure  is  initiated  at  the 
projectile  nose  cone  and  causes  ram  compression  of  the  reactants  via  a  sequence  of  oblique  shocks.  The  shodc 
system  weakens  as  it  mteracts  with  the  expansion  fan  emanating  from  the  projectile  afterbody  shoulder.  Bdiind  the 
projectile,  a  high-speed  wake  structure  is  obtained,  characterind  by  a  wake  reattachment  shock. 

High-IVessure  Compressibility  Effects 

Ram  accelerator  flowfields  are  characterized  by  extremely  high  pressures,  caused  by  rapid  combustion  of  pre¬ 
pressurized  semi-stoichiometric  luel-air  mixtures  under  near-detonative  situations.  Pressures  in  the  neighborhood 
of  KXX)  atmospheres  are  not  unusual.  Under  such  extreme  pressure  conditions,  high  levels  of  conqiressibility  are 
prevalent  and  must  be  accounted  for  during  analysis. 

The  baseline  version  of  CRAFT  was  designed  for  a  multi-component,  thermally  perfect  (but  calorically 
imperfect)  gaseous  mixture  and  the  perfect  gas  equation  of  state  was  assumed.  In  the  RAM  version,  the  perfect 
gas  formulation  was  replaced  by  a  virial  equation  of  sute,  generalized  for  a  multi-component  mixture  as  summarized 
in  Table  V. 
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Table  IV.  Dynamic  Grid  Computations 


EQUATIONS 

r  Qdv  -  f  Qdv  S  •  ?  dsdt  -  4dvdt 

Tenn  I  Tenn  n  Tenn  m 

Term  I:  PRIMARY  CONSERVED  VARIABLES 

Term  H:  FLUX  TERMS 

Term  IH:  SOURCE  TERMS 

V(t):  TIME  VARYING  VOLUME 

METHODOLOGY 

•  TREAT  GRID  AS  PURELY  GEOMETRIC  QUANTITY,  i.e.,  DEFINITION 

OF  CONTRAVARIANT  VELOCITY  UNCHANGED 

•  EVALUATE  FLUX  TERM  BY  ASSUMING  GRID  IS  HELD  CONSTANT  AT 
TIME-AVERAGED  VALUE 

•  INCLUDE  ADDITIONAL  FLUX  ARISING  FROM  CELL  FACE  MOTION 

•  EVALUATE  SOURCE  TERM  FOR  THE  GRID  AT  THE  NEW  TIME 

LEVEL,  THUS  GIVING: 

|q“*‘  -  Q“|^,.,  ^  ♦  F;*‘sAt  *  1 

F,  >  S-P.  and  S  is  the  time  averaged  metric 

NUMERICS 

•  SPECIFY  GRID  MOVEMENT  AS  A  FUNCTION  OF  TIME 

•  EVALUATE  F„  AS  A  THIRD-ORDER  UPWIND  BIASED  FLUX 

•  SOLVE  LINEARIZED  IMPUCTT  OPERATOR  USING  ADI  PROCEDURE 

•  USE  NEWTON  ITERATION  FOR  HIGHER  ORDER  ACCURACY  IN  TIME 

Table  V.  Description  of  Gas  Upgrades 


THE  VIRIAL  EQUATION  OF  STATE  WAS  INCORPORATED  FOR  A  MULTI-COMPONENT  GAS  MDCTURE 


THE  SECOND  AND  THIRD  VIRIAL  COEFFICIENTS  FOR  THE  MDCTURE  ARE  RELATED  TO  THE 
INDIVIDUAL  SPECIES  VIRIAL  COEFFICIENTS, 

IN  GENERAL 


*.(r)  -  E 

HOWEVER,  IF  DATA  IS  UNAVAILABLE  FOR  CROSS  COEITICIENTS  SUCH  AS  B,j,  AN  APPROXIMATE 
RELATIONSHIP  IS  ASSUMED. 


WHERE  Y,  IS  THE  MASS  FRACTION 


Shock  tube  studies  indicate  the  importance  of  compressibility  effects  at  high-pressures.  Figure  S  exhibits  a  10/1 
shock  tube  solution  for  three  different  pressure  ratios  (10/1  atm,  100/10  atm  and  1000/100  atm).  Note  the 
differences  in  the  shock  propagation  pattern  at  the  higher  pressures.  Shock  speeds  are  faster  at  higher  pressures. 
The  steady-state  case  of  Ae  previous  subsection  was  studied  again  with  the  new  RAM  version  of  CRAFT.  Figures 
6  and  7  ^ow  the  con^ruted  contours  of  pressure  and  Mach  number — the  shocks  are  stronger  and  more  closely 
spaced.  Figs.  8  and  9  compare  the  predictions  of  pressure  along  the  tube  wall  and  projectile  afterbody  from  the 
two  computations.  The  effect  of  compressibility  is  evident  and  would  be  more  pronounced  if  combustion  were 
included. 

Adaptive  Grid  Considerations 

The  accurate  prediction  of  transient  ram  accelerator  flowtields  is  strongly  keyed  to  the  ability  to  dynamically 
adapt  the  grids  to  resolve  the  embedded  shock  discontinuities,  wall  boundary  layers  and  separated  ^ear  layer 
regions.  Smart  logic  is  required  to  have  grids  evolve  with  the  transient  strucmre.  Research  has  been  initiated  in 
this  direction  using  solution  adaptive  gridding  based  on  the  NASA  Ames  SAGE  adaptive  solver.^^  Applications  of 
adaptive  gridding  with  transient/unsteady  NS  computations  require  significant  experience  with  regard  to: 

a)  when  to  implement  the  adaptive  solver,  e.g.,  after  how  many  iteration  steps; 

b)  the  associated  adaption  parameter,  e.g.,  a  certain  weighted  combination  of  pressure,  Mach 
number,  temperature,  etc.;  and, 

c)  the  control  parameters,  e.g.,  the  degree  of  grid  distortion,  orthogonality,  etc.,  per  adaptive  pass, 
etc. 


The  development  of  a  suitable  grid  adaption  strategy  for  ram  accelerators  has  been  initiated  with  focus  directed 
towards  developing  rules  for  some  of  the  considerations  outlined  above.  Preliminary  adaption  strategies  are  being 
developed  for  a  steady-state  trajectory  point  aiulysis  and  will  be  subsequently  implemented  for  transient  calculatioiu 
in  a  *hands-ofr  fashion.  Figures  10  and  11  display  examples  of  adapted  grids  generated  for  the  steady-state  case 
analyzed  earlier.  The  grids  were  generated  by  selecting  varied  a^ption  parameters  and  controlling  the  grid 
distortion.  Figure  12  shows  the  grid  and  resultant  solution  after  several  passes  of  the  solution  adaptive  grid 
generator  and  fortber  solution  convergence  with  the  updated  grids.  Compare  the  crispness  of  the  shock  waves  with 
the  baseline  grid  predictions  of  Figs.  3  and  4. 

Finite-Rate  Chemistry  Considerations 

Finite-rate  chemistry  mechanisms  and  rates  are  required  for  CFD  simulations.  Unforturutely,  the  vast  rrujority 
of  hydrocarbon  kinetics  data^^’^^  were  acquired  (and  mechanisms  formulated)  under  laboratory  conditions  with 
limited  (if  any)  applicability  to  the  high  pressure,  semi-detonative  ram  accelerator  environment.  Extrapolation  of 
kinetics  data  to  higher  pressures  may  lead  to  numerical  predictions  of  dubious  quality.  Under  such  constraints, 
global  multi-step  mechanisms  provide  a  useful  compromise  until  more  elaborate  mechanisms  can  be  formulated  via 
the  use  of  computational  chemistry,  as  was  demonstrated  recently  for  NASP  in  the  hypersonic  scramjet  environment 
by  Cooper  and  coworkers  at  NASA  Ames. 

One-dimensional  shock  tube  studies  at  laboratory  conditions  have  been  initiated  to  gain  a  prelimirury 
understanding  of  the  kinetic  processes  relevant  to  the  ram  accelerator.  The  objective  of  these  simulations  is  to 
identify  the  interaction  between  shock  formation  and  ignition  delay  in  a  transient  environment.  A  related  objective 
is  to  also  establish  an  upper  pressure  bound  for  the  applicability  of  conventional  hydrocarbon  kinetics  rates  and 
mechanisms.  Figure  13  shows  the  schematic  of  a  problem  studied  and  the  kinetic  mechanism  utilized  for  the 
simulations.  Figures  14  and  IS  show  the  temperature  and  pressure  distribution  in  the  shock  tube  at  two  instances, 
separated  by  S/is.  At  the  earlier  time,  the  shock  front  is  clearly  seen  (characterized  by  the  von  Neumarm  spike), 
followed  by  a  constant  pressure  induction  zone  and  subsequently  the  flame  with  rapid  temperature  rise.  The  flame 
front  travels  faster  than  the  shock  and  merges  with  it  a  little  later,  as  shown  by  the  relatively  uniform  post-shock 
temperature.  Figures  16  and  17  show  the  distribution  of  major  and  minor  species  at  the  two  time  instances 
considered  and  a  consistent  behavior  is  observed.  Further  one-dimensional  studies  with  upgrades  rates/mechanisms 
are  felt  to  be  essential  for  gaining  an  understanding  of  what  might  occur  under  actual  ram  conditions. 
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CONCLUDING  REMARKS 


This  paper  has  described  the  CRAFT/RAM  computer  code  developed  for  application  to  ram  accelerator 
flowfields.  The  code  developed  is  a  specialized  version  of  the  CRAFT  research  code  with  generalized 
thermochemistry,  advanced  turbulence  modeling  and  specialized  dynamic/solution  adaptive  gridding  capabilities. 
Issues  pertinent  to  developing  a  fundamental  understanding  of  the  detailed  physical  and  thermochemical  processes 
governing  ram  accelerator  concepts  were  exhibited  with  numerical  simulations  illustrating  the  code’s  capabilities 
to  simulate  these  processes.  Future  work  will  include  boundary  layer  effects  (projectile  and  launcher  tube),  full 
transient  simulations,  and  exploration  of  thermochemical  issues. 
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Figure  2.  a)  Baseline  301  x  71  curvilinear  body>filted  grid,  and  b)  Expanded  view  of  grid 
baseline  region. 


Figure  3.  Computed  pressure  contmirs  for  SSimn  projectile. 


Figure  4.  Computed  Mach  number  contoun  for  38mm  projectile. 
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Figure  9.  Con^Mrison  of  projectile  surftce  pnesure  distiibutioa 
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12a. 
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Figure  12.  «)  Grid  obtaiimd  after  multiple  adaptiooa;  b)  Pfcsaare  ooutoura  on  adqMed  grid 
aad  beaeline  grid;  and  c)  Madi  ouinber  cootoun  oo  ^apled  grid  and  beaeliae  grid. 
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ONE-DIMENSIONAL  SHOCK  TUBE 


METHANE  REACTION  MECHANISM 


Reaction 


kf,  rate  constant 


1.  CH4  +  M-CH3  H  -t-  M 

2.  CH4  +  Oj  -  CH3  +  HOj 

3.  O2  M  -  20  M 

4.  CH4  +  0  -  CH3  +  OH 

5.  CH4  +  H  -  CH3  +  Hj 

6.  CH4  +  OH  -  CH3  +  HjO 

7.  CH3  +  0  -  HjCO  +  H 

8.  CH3  +  O2  -  H2CO  +  OH 

9.  H2CO  OH  -  HCO  +  H2O 

10.  HCO  -t-  OH  -  CO  +  H2O 

11.  CO  +  OH -002  +  H 

12.  H  +  O2-OH  +  0 

13.  0  +  Hj-OH  +  H 

14.  0  <«•  H2O  -  20H 

15.  H  H2O-H2  *  OH 

16.  H  +  OH  +  M-H2O  +  M 

17.  CH3  +  O2-HCO  +  H2O 

18.  HCO  -t-M-H  +  CO-^M 


1.5  X  10’®  exp(-100,600//?7) 

1.0  X  10’*  exp(-45.400/R71 

3.6  X  10’“r  ’  ®exp(-118.800//?r} 

1.7  X  10’®exp(-8760/R71 

6.3  X  10’*  exp»-12.700/R71 

2.8  X  10’®  exp(-5000/R7) 

10’®  -  10’® 

10”  -  10’* 

10’®  -  10’® 

10’2.  10’® 

3.1  X  10”  exp(-600/R7) 

2.2  X  10’*  exp(-16,600//?7) 

4.0  X  10’*  exp|-9460/R71 

8.4  X  10’*  expi-18.240//?n 
1.0  X  10’*  exp<-20.400//?n 
2.0 X  lO’r*’® 

to”  - 10’® 

2.0  X  10’®r’'2exp(-28,800//?71 


Figure  13.  Reactive  shock  tube  schematic  and  methanr.  Idnetic  mechanism 
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Figure  14.  Temperatun  distributioa  in  shock  tube 
at  two  tune  inslanoee. 
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Figure  15.  Pleasure  distributioa  in  shock  tube 
at  two  time  instanoee. 
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Figura  K.  Major  q>ecies  distributiaA  io  diode  tube  at  two  tune  instancea. 
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Figure  17.  Minor  speciea  diatributioo  in  diock  tube  at  two  time  inatancea. 
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